Weakly localised bosons 



D. Das and S. Doniach 
Department of Applied Physics 
Stanford University, Stanford, CA 94305. U.S.A. 
(Revised Draft 21 Jan 98) 

Abstract 

In this paper, we argue that, in 2-d, the weak localisation of bosons, which 
occurs on the insulating side of the superconductor-insulator transition, is 
characterised by p ~ Zn(l/T), as compared to a ~ InT for fermions. Such 
an unconventional behaviour is tied to the diffusion pole of the delocalised 
(uncondensed) vortices. 
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Superconductivity in low- dimensional systems has been of considerable interest over the 
past few decades. This has received a strong boost with the discovery of high-T c superconduc- 
tivity. High-T c superconductors are layered materials and exhibit a superconductor-insulator 
transition at low doping. In an attempt to understand the two-dimensional aspect of this 
new phenomenon, an appreciable amount of effort has been focussed on the low-temperature 
behaviour of thin films made out of the low-T c superconductors [1]. In these experiments, 
one observes that the film undergoes a transition from a superconductor to an insulator at 
T — > 0, as one tunes the thickness, disorder or the magnetic fieldfl]. There are two dis- 
tinctively different schools of thought, for explaining this set of observations. One school 
[2] believes that the insulating behaviour occurs due to the localisation of charged bosons, 
viz. Cooper pairs. Whereas, the other school [3] believes that during the said transition, 
Cooper pairs break up into fermions which naturally localise in 2-d and cause insulating 
behaviour. Each school is supported by a corresponding theory. Cooper pair localisation 
is supported by Fisher's theory[4], whereas the fermion localisation picture is interpreted 
in terms of Fukuyama's theory [5]. In the latter scenario, Cooper pairing is hampered as 
a result of suppression of density of states due to enhanced coulomb repulsion. Although 
Fisher's theory seemed to work in many cases [6], there has been unflinching evidence[3,7] 
that Fukuyama's theory is right for certain systems. Thus each mechanism probably has an 
appropriate class of materials for which it is relevant. 

So far, there has been no straightforward way to find which theory applies in a given 
system. In this paper, we suggest a simple way to resolve this debate. In particular, we 
argue that, in 2-d, weak localisation of bosons implies p ~ ln(l/T), whereas for fermion 
localization one expects a ~ InT [8]. The closest experimental observation of this non- 
fermionic effect seems to be the experiment by Gerber on granular Pb films [9]. A similar 
behaviour has been reported in Ref. [6]. It is interesting to note that such a p ~ ln(l/T) 
dependence has been observed in the high-T c materials[10] as well, where bosonic models 
have been suggested to be relevant [11]. In what follows, most of our discussion will be 
restricted to two dimensional systems at finite temperatures, the focus being exclusively on 
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low-T c superconductors. 

Why not a ~ InT for bosons? a ~ InT behaviour is profuse in 2-d electronic systems. 
Many of these logarithmic tendencies are associated with the critical properties of electrons, 
viz. the lower critical dimension of the electronic system is two. By contrast, in a bosonic 
system, the lower critical dimension is one [4]. There is a further difference between the 
two kinds of charge carriers. Let us consider, for example, a system of non-interacting 
electrons. Here a ~ InT follows from one-parameter scaling[8]. It has been argued that 
[12] the scaled conductance g = G/(e 2 /h) for such systems depends on a single parameter 
AE/5W, where AE = bandwidth and 5W = (iVo-L^) -1 is the scatter in the energy levels 
in a system of size L in d-dimensions. N is the density of states at the Fermi level. Let 
us repeat the argument [12] here for the sake of completeness. In this picture, one considers 
constructing a sample of size (2L) d , in d dimensions, out of 2 d samples of size L d . A given 
eigenstate of the larger sample is a linear combination of the eigenstates of the smaller 
samples. What really matters here is the amount of admixture among the (old) eigenstates, 
which in turn is controlled by the overlap integral and energy denominator. The measure 
of the energy denominator is 8W . To obtain an estimate for the overlap, one considers 
repeating a given L d sample in one direction, subject to appropriate boundary conditions. 
Each eigenstate broadens to form a band, and the bandwidth AE is an estimate of the 
overlap integral. If the wave function is localised, AE, and hence, AE/5W, is exponentially 
small. If the wavefunction is delocalised, AE/5W is large and AE is sensitive to boundary 
conditions. Thus, the scaled conductance g is a simple function of a single parameter 
AE/5W. But this is not enough for a bosonic system. The phase of the wavefunction is 
an important parameter here. A pretty argument due to Anderson[13] suggests that, for 
the Bose case, there are strong phase relationships among the (smaller) samples, discussed 
above, which would invalidate one-parameter scaling and also the logarithmic behaviour of 
the conductivity along with it. More precisely, in the localised phase, the so called Bose- 
glass(BG) phase[14], there are big patches of locally superfluid regions which strongly modify 
the scattering of the charge bosons in the following way. When a charge boson scatters off 
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an impurity, it correlates with(or, talks to) all the charge bosons within a distance of the 
order of the localisation length £; oc . As a result, the boson scatters as a blob, i.e. the 
scattering is inherently multiparticle/collective in nature, as compared to a single particle 
scattering typical of a non-interacting particle. (Plus, this is also the physical reason why 
we find a tenfold enhancement of resistance over the fermionic counterpart.) That is why it 
is much more convenient to think in terms of the collective excitations, viz. vortices, which 
are dual to the charge (boson) degrees of freedom [15]. Translated into the vortex picture, 
big superconducting patches convert into regions where vortices are localised. In between 
these large regions, there are narrow channels where the vortices are delocalised, and only 
these make a significant contribution to the vortex conductivity [16]. These channels are 
superfluid below a certain temperature T AV (please see later), above which they lose phase 
coherence [17]. Further, we note that [18] the charge resistivity p c is related to the vortex 
conductivity a v by 

p c = (h/4e 2 )a v (1) 

In the rest of the paper, we are going to exploit this duality relation between charges 
and vortices extensively and evaluate the resistance as seen by the charges in terms of the 
vortex mobility. 

In order to calculate the vortex conductivity, or equivalently the charge resistivity, we 
consider a model of interacting (charged) bosons in a random potential[14]. Before we write 
down this model explicitly, let us discuss the physical picture [Fig. 1] generated by this model. 
Consider a system of bosons(e.g. Cooper pairs in a granular superconductor) interacting 
with a short-range repulsive potential V in the presence of a finite disorder A. Let J be 
the zero point energy of the bosons. When the ratio (J/V) is large, the bosons condense 
into a superfluid(SF) state. As the ratio (J/V) is decreased (or, alternatively, (A/V) is 
increased), the bosons tend to localise and superfluidity decreases. At a critical value of the 
kinetic energy (J/V) c , (global) superfluidity is completely destroyed and the bosons make 
a transition to the Bose-glass phase, where the bosons are essentially localised by disorder. 
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On the insulating side, very close to the phase boundary, the bosons are weakly localised 
and there are finite regions (of the order of the localisation length £ Joc ) where the bosons 
form a local condensate. As one moves further from the phase boundary by reducing (J/V), 
one gradually enters the strongly localised regime. Here the overlap between the nearest 
neighbours is really small and Mott variable range hopping is expected to occur [18]. In this 
paper, we address the issue of conductivity of Cooper pairs in the Bose glass phase very 
close to the BG-SF phase boundary. For the sake of completeness, we note that when the 
average number of bosons per site is an integer, there is an additional Bose insulating(BI) 
phase. In this phase, there are exactly n bosons per site (where n is an integer) and the 
phase is gapped. 

One can consider this entire phase diagram in the (dual) vortex picture as well. In the 
superfluid region, the vortices form closed loops in zero field. At the superfluid-Bose glass 
transition, the loops blow up as a result of the quantum zero point motion at T = and 
thermal motion at finite T, and the vortex-antivortex pairs break up [19]. Since these vor- 
tices and antivortices are bosons, they form a bose condensate at T = 0[15]. However, for 
T > T AV , the lambda transition of the superfluid vortex condensate [20], this is destroyed, 
because the disordering of the charges(viz. Cooper pairs for this case) [18], which represent 
the topological excitation over the vortex degrees of freedom, will destroy any possible phase 
coherence among the vortices. We further note that the existence of a vortex condensate 
would simply mean, from (1) that, p c = oo at a finite temperature. This signature of a 
"quantum insulator" is rather counterintuitive, so far as the resistive behaviour of classical 
insulators is concerned. However, it can be made clearer in the dual "cooper lattice" phase 
where the system is essentially a charge density wave in which the pairs are strongly cor- 
related. The superinsulator behaviour then corresponds to the response of a sliding charge 
density wave in the absence of pinning. This behaviour would probably be suppressed ex- 
perimentally by finite size and nonlinear effects. In the rest of the paper, we shall work in 
the regime T > T AV which we assume may be quite small in the materials such as those 
studied by Gerber. 
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Although the above phase diagram is simplest to model on a lattice, in the rest of this pa- 
per we shall consider the bosons in a continuum [2 1,1 5, 2 2]. The main purpose of this exercise 
is to demonstrate that under the duality transformation, vortices see potential disorder(as 
seen by the charges) as a random magnetic field(RMF). This feature is particularly explicit 
in a continuum description. Hence we consider the following model Hamiltonian(/i = 1) 

H = J d 2 x[^ | ±fy | 2 +\ J d 2 y5p{x)V{x - y)5p{y) - p(x)p(x, r)], (2) 



where 5p = p — p, p being the (neutralising) background charge density with partition 
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function Z = J Vi/j*Vi/jexp(-S), and S = J drd 2 xC = J dr[J d 2 xip*^ + H}. Here V(x) is a 



short-range potential V(x) = VS(x). 

Equivalence with the above phase diagram is established by setting p/m = J. Now, 
invoke a duality transformation [23] ip = ^/pe l6 (j) vi where 6 = non-topological part of the 
phase and (f> v =topological part of the wave-function, with = 1. Also, p = po + 5pi(x), 
where Spi(x) represents the impurity potential which is Gaussian distributed with variance 
A. (We have not written down the chemical potential po explicitly in what follows for 
convenience, and have taken p(x) as 5pi(x).) Further, define a gauge field Mi(x) (i = 1,2) 
by the relation e^diMj = —p{x)/V = B(x). Here M(x) is independent of r. As a result, 
one is led to the following Lagrangian - 

Tfi —* ~ — 

C = 2p 1 J |2 + 2" (5p)2 + 2mj ^ + 2m ^ M k, (3) 
(p = 0, l,2,k= 1,2), with 

< B(x)B(x>) >= (A/V) 2 5(x-x'). 

Here J M = (Sp, J) = e^ vp d v a p represents the charge (boson) density and current, whereas 
iii = j p = (1/2^)6^^(0*19^01,) represents the vortex density and current[24]. Integrating 
out the charge degrees of freedom, and transforming over to the real time, we are led to the 
following model Hamiltonian 

H v = ^2 7T—(Pia ~ 2nqiM ia ) 2 + \^n 2 — ^ ft&MI x i ~ x i I /Co) + consts (4) 



Here % = ±1, (a = 1, 2) represents the charge on a vortex and an antivortex respectively, X{ 
and Pi represent the positions and momenta of vortices and summation is over all vortices 
and antivortices. Here m v refers to the vortex mass. Details of the notation and missing steps 
can be found in [23]. Thus, we find that vortices move in a random effective magnetic field 
under a duality transformation. A few comments are in order. In going from eqn (3) to (4), 
we have included a vortex mass term. This kinetic energy term is essential for the quantum 
melting of the vortices and is guaranteed by the underlying electronic degrees of freedom[25]. 
Also, we have neglected a current- current interaction term among vortices resulting from the 
short-range interaction among charges. This is valid when n v /Vm v <C l/47i 2 (n v =vortex 
density) and at low temperatures. The physical content of this condition is as follows. The 
mass term for vortices generates a long-range interaction among the charges[26]. The above 
inequality simply states that this long-range component is much weaker than the short-range 
component. Thus, the BG phase discussed above is perturbed to as little extent as possible. 

As we noted earlier, in the BG phase, the vortex-antivortex pairs break up. Since they 
are delocalised, they screen each other, and as a result, interact via short-range interactions. 
Thus, we are led to evaluate the conductivity of a vortex liquid as it diffuses in the presence 
of a random magnetic field at finite temperature. We shall do this using perturbation 
theory. As is well known in the context of the fermionic problem[27,28], the diffusion pole, 
viz.r(g,cj„) = (l/mr 2 )/(| u n \ +Dq 2 ), for e m (e m + uj n ) < 0, (where u n = 2iinT, n =integer, 
are Matsubara frequencies and D is the diffusion constant) generates singular corrections 
to the conductivity(at finite temperature). In the absence of a condensate, this is true of 
the Bose case as well, with m = m v ,D = jv 2 r tr , where \m v v 2 — /i v = =chemical 
potential of vortices, r tr =transport time in a RMF= m v £ 2 /iT 2 (A/V) 2 , and £ =pair size, 
defines the microscopic scale of the Bose problem, r = the elastic scattering time in the 
RMF, and we will take r = r tr [28]. Now, for superconductors, fi v = d(<& /4:iT\) 2 lnK, where 
k is the Ginzburg-Landau parameter and d is the film thickness. One can rewrite this 
as ji v = aop/m = aoJ, with a = jinn. Using the usual value for the Josephson coupling 
energy J = (i?g/2i? n ) A , where A is the superconducting gap, Rq = h/Ae 2 ~ 6.45.KX2, and 
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R n is the normal resistance of the sample, we obtain the condition for validity of perturbation 
theory as 

2 k F d E F Rq 

Because magnetic disorder is time-reversal symmetry breaking, the cooperon mode is sup- 
pressed^]. The diagrams which contribute are the Altshuler-Aronov type diagrams, shown 
in fig. 2. We merely quote the result here and refer the reader to Refs. [27,28] for more 
details. We obtain 

8a v = -A ln(Tr t r) (6) 

an extra minus sign coming from the replacement of anticommutators by commutators 
for bosons. Physically, this makes sense because as we go down in temperature the vor- 
tices being bosons tend to bose-condense. Here A = (2/n)(l + F/2), a number of order 



unity, with F = \/yJ\ + (a /27r 2 ). Eqn.(6) is an enhancement over a background part 
Rvo \rmf= (h/4e 2 )(nT tr /m v ) ~ {h/Ae 2 ){n v ^l/ {A/V) 2 ). The constant part of the conductiv- 
ity will receive contributions from the Bardeen-Stephen processes, R B s = 2vri? n (n 1 ,^o)[29]. 
Putting all these contributions together, we obtain, 

Pc = R b + Roln{T OB /T), (7) 

where R b = [2nR n + (R Q /(A/V) 2 )](n v Q), and i? = A R Q . Also, in the BG phase, we 
expect (A/V) 2 to be a number of order unity. Here T 0B = l/r tr . Eqn.(7) is the main result 
of the paper. 

Comparison with experiment. The logarithmic behaviour for bosons, as predicted by 
eqns.(6) and (7), is much steeper than the analogous fermionic behaviour (with interactions 
included), viz. a F = a - (2 - 2ln2)(e 2 /nh)ln(T 0F /T) [12], with a = (e 2 /h)(k F l). For fairly 
large values of k F l, this implies 

PF ~ p + R F ln(T 0F /T), (8) 
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where 1/Tqf = tf =electronic transport time, and Rf = (2 — 2ln2)(p /nkFl)- Since we do 
not have all the necessary material parameters available for Gerber's experiment [9], we make 
some estimates based on Mo-Ge samples discussed in Ref.[6]. In this case, Tav ~ 10mK[20], 
much lower than the temperature range of measurements. This justifies our neglect of the 
"quantum insulator" phase in this discussion. Taking k F l ~ 5 [30], we have R F ~ 0.1 — IKVl. 
The scale of Ro is set by Rq = (h/4e 2 ) and we have Ro ~ 5kQ. Thus, we expect a 
change in slope of a R vs InT curve by a factor of about 10 when bosonic conduction sets 
in. There is one more difference between Bose and Fermi resistivities. In the case of the 
former, Tqb — 0.1 — 10i^, whereas Tqf — 10 2 — 10 3 i^. Thus, the fermionic mechanism is 
a high temperature phenomenon, whereas the bosonic behaviour is a predominantly low 
temperature phenomenon. For low-T c materials, T 0B ~ T p , the pair formation temperature, 
and the InT dependence may be limited by T p and set in as soon as the pairs are formed. 
These features are clearly observable in Gerber's experiment [9] . Also, in this experiment, at 
a very high magnetic field which kills quasi- reentrance (characteristic of pair formation), this 
low temperature logarithmic behaviour disappears as well. We consider this to be a telltale 
evidence of bosonic transport. However, it must be mentioned that in this experiment [9], 
the scale of the resistivity is unusually high, lying between megaohms and gigaohms. For 
a typical resistivity of fig. 2 in Ref.[9], we get Rq/R u ~ 10~ 4 . This means, from eqn(5), 
that E b r tr ~ 0(1). So, higher order corrections in (1/E b r tr ) are particularly important 
in this case and probably a certain class of diagrams needs to be resummed. We also 
take this opportunity to comment that Wolf's observations [31], which originate on the high 
temperature side, are most likely not due to bosonic transport, but have some fermionic 
mechanism attached to it. 

Our focus in this paper has been only on short range interactions among the charge 
bosons, although in the superconducting films like Mo-Ge[6], which are of current interest, 
long-range interactions dominate. It is quite possible that inclusion of such effects will 
renormalise the coefficients only, rather than affecting the temperature dependence strongly. 
This needs further investigation. 
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All these issues prevent us from making good contact with the existing experiments and 
we only hope to have conveyed to the reader the essence of bosonic transport, and how it 
is distinguished from its fermionic analog. We believe, however, that this behaviour should 
be seen in other systems where bosonic transport dominates. 

What this is not? We would like to comment here, before closing, that this model is 
not the same as the traditional Bose Hubbard model with simple on-site repulsion. The 
latter does not generate a vortex mass term [32] which is particularly important for vortex 
mobility calculations discussed above. Thus, any simulation which tries to track this resistive 
behaviour must include an appropriate vortex mass term along with the usual Bose Hubbard 
model terms. 

To summarise, we have argued that in the regime where Cooper pairs are weakly localised, 
the vortices are in a liquid phase. Thermal diffusion of this quantum liquid in the presence 
of (pseudo-magnetic) disorder leads to logarithmic temperature dependence of the (charge) 
resistivity. This kind of temperature dependence should be observable in a bosonic system, 
so long as the temperature is higher than that at which the vortices themselves would form 
a bosonic condensate. 
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FIGURES 

FIG. 1. A schematic phase diagram for the Bose localisation problem. The hatched region 
shows the region of interest. 

FIG. 2. Altshuler-Aronov diagrams 
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Fig. 2 (BLR685) 




